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We investigate the local density of states, spectral function, and superconducting pair amplitudes for 
signatures of Majorana fermions in vortex cores in ferromagnetic and spin-orbit coupled semiconductor- 
superconductor heterostructures. We show that the Majorana fermion quasiparticle momentum distribution is 
always symmetrically distributed at a finite radius around a high symmetry point, thereby providing a necessary 
condition for a low-energy state to be a Majorana fermion. In real space profiles of the local density of states 
through the vortex core the Majorana fermion, together with other finite-energy vortex states, form a characteris¬ 
tic x-shape structure only present at non-trivial topology. Moreover, we find that the Mexican hat band structure 
property of the topologically non-trivial phase translates into multiple high-intensity band edges and also vortex 
core states located above the superconducting gap in the local density of states. Finally, we find no strong cor¬ 
relation between odd-frequency pairing and the appearance of Majorana fermions, but odd-frequency pairing 
exists as soon as ferromagnetism is present. In fact, we find that the only vortex superconducting pair amplitude 
directly related to any phase transition, is the appearance of certain spin-triplet p-wave pairing components in 
the vortex core at a pre-topological vortex core widening transition. 

PACS numbers: 74.90.+n, 03.65.Vf, 74.20.Rp, 74.25.Uv, 74.55.+V 


I. INTRODUCTION 

In recent years the study of topological phases has lead to 
the prediction that so-called Majorana fermion quasiparticles 
can appear in certain types of topological superconductors .G® 
These Majorana fermions are of great interest for two quite 
distinct reasons. First of all, Majorana fermions are hypo¬ 
thetical particles long sought in particle physics,but so far 
without any conclusive evidence in favor of their existence.!^ 
Majorana fermion quasiparticles in the solid state, although 
not fundamental particles, are analogous to their fundamental 
counter-parts, and may therefore provide an independent way 
to discover Majorana fermions!^ The second reason they are 
of interest is because when they emerge in condensed matter 
systems, typically localized on various defects, they do so in 
a way which makes the ground state degenerate. These dif¬ 
ferent ground states can in two dimensions be continuously 
deformed into each other by braiding the defects around each 
other. Such operations are predicted to be non-Abelian, and 
braiding Majorana fermions may thus be utilized for the im¬ 
plementation of robust topological quantum computing.® 

One place where Majorana fermions are expected to emerge 
is in vortex cores in certain two-dimensional (2D) topolog¬ 
ical superconductors. A prominent example of such a sys¬ 
tem is provided by a heterostructure of a thin layer of Rashba 
spin-orbit coupled semiconductor, sandwiched between a fer- 
romagnet and a conventional s-wave superconductor.®!! Ma¬ 
jorana fermions are also expected t o occ ur at the end points 
of ID wires of similar composition.E^E^ Possible signatures 
of Majorana fermions has already been reported f or suc h ID 
systems, n^^although some results are still debated.^^EHlEven 
if braiding in principle is possible for ID systems through the 
use of wire networks,!^ vortices lend themselves more natu¬ 
rally to be braide d, fo r example through the use of magnetic 
force microscopy.!^®! In the light of the potential versatility 
of vortex Majorana fermions it is of large interest to make a 
thorough investigation of different types of signatures for Ma¬ 


jorana fermions in superconducting vortices. 

Analytical and numerical results on vortices have al¬ 
ready predicted that Majorana fermions appear in the 
topologically no n-triv ial phase of the above mentioned 
heterostructures in this work we carefully investi¬ 

gate experimentally relevant signatures of both the Majorana 
fermions and the topologically non-trivial phase. In particu¬ 
lar, we focus on signatures in the spectral function, local den¬ 
sity of states (LDOS), and superconducting pair amplitude. A 
simultaneous investigation of these three quantities is benefi¬ 
cial, as many signatures in one or another of these are closely 
related to particular features also in the other properties. 

More specifically, we find that the Bogoliubov-de Gennes 
(BdG) quasiparticle spectrum for a Majorana fermion is 
strictly required to be symmetrically distributed at a finite ra¬ 
dius around a high symmetry point, which thus provides a 
necessary condition for any candidate Majorana state. We 
also find that the vortex Majorana fermion and the generic 
finite-energy Caroli-Matricon-de Gennes vortex stateJ^ are 
well separated in energy. Together they form a characteris¬ 
tic x-shape structure in subgap LDOS profiles right through 
the vortex core only in the topological phase. The Majorana 
mode is well-localized to the center of the core, while the 
finite-energy states disperse further out from the center. Be¬ 
yond the occurrence of Majorana fermions, the topologically 
non-trivial phase can also clearly be distinguished by a Mexi¬ 
can hat shaped band structure, which can be probed with mo¬ 
mentum space probes such as angle-resolved photoemission 
spectroscopy (ARPES). The Mexican hat shaped band struc¬ 
ture further gives rise to multiple band edges, showing up as 
double peaks in the DOS. The existence of the double peaks 
enables experimental techniques sensitive to the DOS to act 
as a probe of the topological phase, possibly most relevant for 
scanning tunneling spectroscopy (STS). The two band edges 
in the double peak also behave differently in the presence of 
a vortex. While one edge collapses to give rise to a Majorana 
fermion and the generic Caroli-Matricon-de Gennes subgap 
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vortex core states, the other band edge is instead pushed up in 
energy and there gives rise to a second set of vortex core states 
that appears as a rising band edge in the LDOS spectrum of the 
vortex core. We note that these types of STS measurable sig¬ 
natures are of particular interest, as the feasibility of such ex¬ 
periments recently have been demonstrated in a related setup 
consisting of an s-wave superconductor with vortices coated 
by a topological insulator.lSl Finally, we show that there is no 
distinct onset of odd-frequency pairing as a result of the ap¬ 
pearance of Majorana fermions, but odd-frequency pairing is 
present as soon as there is finite magnetism, independent on 
the topological phase. However, we find a strong correlation 
between a pre-topological vortex core widening transition and 
the onset of p-wave pair amplitudes in and around the vortex 
core. 


II. MODEL 

We here consider a 2D topological superconductor with 
the essential building blocks being s-wave superconductiv¬ 
ity, Rashba spin-orbit interaction, and Zeeman ferromagnetic 
term. To achieve self-consistent microscopic details for a vor¬ 
tex core we study the system on a square lattice. The free pa¬ 
rameters of the model are then the nearest neighbor hopping 
t (setting the kinetic energy), chemical potential p, Zeeman 
field Vz, Rashba spin-orbit interaction a, and superconduct¬ 
ing pair potential 14c- All energies can be measured relative 
to the kinetic term, which we do by setting t = 1. The H amil- 
tonian describing this system can be written 

H = Hkin+'Hv^ +'Hso +'Hsc, ( 1 ) 

'Hkin = ~t 'y ' ~ M ^ ' Cj^CicT, 

(ij),cr 

'^Vz ^ ^ aa'■> 

'^SO = - f 

i 

+i(4_y;Cit - -f H.c.J , 

Tisc = ^ H.C.). 

i 

Here i and j are site indices on the square lattice, a is the 
spin index, and c|^^(cicr) is the electronic creation (annihila¬ 
tion) operator. We are primarily interested in a lightly hole- 
doped semiconductor, which is achieved by setting p, = 4. 
The superconducting order parameter Ai enters as a param¬ 
eter in the Hamiltonian, but is determined self-consistently 
using a superconducting pair-potential 14c- This is done by 
solving Eq. Q within the BdG formalism, and re-calculating 
the order parameter using 

a'”*'' = - 

= -r.. E ( 2 ) 

E ^<0 



FIG. 1: (Color online). Energy spectrum for a superconducting vor¬ 
tex as a function of the Zeeman field (red and green lines), strength of 
the superconducting order parameter (thick blue line), and absolute 
value of the Zeeman field (black dashed line). The I, F, and II regions 
corresponds to the trivial phase, trivial phase with wide vortex core, 
and topologically non-trivial phase, respectively. Four representative 
sample points are also marked, used in other figures. 


Here is the electron up (hole down) component on 

site i, and m is the iteration step. We are able to study a single 
vortex by specifying the initial order parameter configuration 
0) = |A(°)(r, 0)|e“*®, but then letting the supercon¬ 
ducting amplitude and phase fully relax. This allows for a 
fully self-consistent order parameter profile to be obtained un¬ 
der the single requirement that the phase winds 27r around the 
vortex core. 


A. Phase transitions 

We have previously shown that close to the topo¬ 
logical phase an unrelated phase transition can take 
place, due to the competition between ferromagnetism and 
superconductivity.!^ Xhis phase transition manifests itself in 
the sudden widening of the vortex core, together with a jump 
in the magnetization in the core. This means that the phase 
diagram for the system can be divided into three different 
regions which we label I, T, and II. These are the topologi¬ 
cally trivial, topologically trivial but with a wide vortex core, 
and topologically non-trivial phases, respectively. In Fig. [T] 
the energy spectrum as a function of the Zeeman ferromag¬ 
netic field is reproduced, and the three phases are marked with 
corresponding labels. In the topologically non-trivial phase 
II, zero-energy Majorana states appear (red line). However, 
states at or close to E = 0, which are not Majorana fermions. 
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also appears in the F region. It is therefore directly clear that, 
if looking for a single signature such as a state at i? = 0, there 
is a significant risk of mistakenly identifying a state as a Ma- 
jorana fermion even though it is not. In this work we therefore 
carefully investigate several different experimental signatures 
of both the Majorana fermions directly, as well as signatures 
related to the different phases. Most of these signatures can 
be directly accessed with real space and band structure probes 
such as STS and ARPES, respectively, while a few other sig¬ 
natures are of at least important conceptual value. 


the topological phases of the system. First of all, because 
we study a lightly hole-doped semiconductor, the relevant 
condition for being in the topologically non-trivial phase iP 
(4f — -f |Ap < 14 ^ < fj? + |Ap. We note that in the 
bulk there is no difference between the two topologically triv¬ 
ial phases I and F. The difference between these two phases 
only become apparent in self-consistent vortex calculations. 
Further, the band structure of the system is given bjP 


III. BULK BAND STRUCTURE 


Before turning to the results for a vortex, we begin with 
a few important remarks about the bulk band structure and 


En (k) = ± ye2(k) + a 2 /: 2 (k) + + |a|2 ± 2 ^e 2 (k)a 2 £ 2 (k) + (e2(k) + |A|2) 442 , (3) 


where n is the band index, e(k) = —2t {cos{kx) +cos{ky)) — 
fji the kinetic energy, and £o(k) = {sm(ky), — sm{kx)) the 
spin-orbit coupling. In Fig. we plot the band structure for 
representative points of the trivial and non-trivial phases, as 
well as points where the band gap closes. It is clear from the 
plot for phase II that in the topologically non-trivial phase the 
two bands closest to E = 0 form two oppositely facing Mex¬ 
ican hats centered around the point k = (tt, tt). [More gener¬ 
ally this could also be at the other high symmetry point; e.g. in 
a lightly electron-doped semiconductor the Mexican hats are 
centered around k = (0,0)]. The trivial band structure in 
phase I instead takes the form of two ordinary parabolas. Fur¬ 
thermore, the two plots (D) and (M) shows two different ways 
through which the bulk band gap of the non-trivial phase can 
be closed. At the topological phase transition (I •(-> II), the 
band gap closes by forming a Dirac cone. On the other hand, 
the band gap can also be closed by letting A —> 0, in which 
case the system becomes metallic. In the top left figure a 
schematic view of the two different routes through which the 
gap can be closed is shown. From a topological point of view, 
the only important thing is that the I and II regions are sepa¬ 
rated by a gap closing, and the two types of gap closings can 
be considered topologically equivalent. We will however see 
that by distinguishing between the two types of gap closings, 
the origin of important experimental signatures of Majorana 
fermions as well as the topological phase can be understood. 


A. Edge modes 

In the semi-classical limit it is possible to treat real and re¬ 
ciprocal coordinates as independent of each other. The edge 
states that appear at interfaces of topologically non-trivial 


phases can then be understood as a consequence of the bulk 
band structure being required to go through a gap closing 
when passing from the topologically non-trivial phase inside 
to the trivial phase outside!^ Zero-energy modes on impuri¬ 
ties can be understood to appear for similar reasons. However, 
in reality interface and impurities are fairly abrupt in nature, 
and it is therefore not obvious that a semi-classical treatment 
relying on the bulk band structure reliably predicts the proper¬ 
ties of such defects. Nevertheless, we know that the prediction 
of zero-energy states at interfaces is still valid also for abrupt 
defects, when the forced gap closing is associated with a jump 
in a topological invariant, such as the Chern number. This 
have been demonstrated for many systems both numerically 
and analytically, and can in some cases be formally justified 
through the use of an appropriate index-theorem, see e.g. M 


1. Semi-classical prediction 

In what follows we will repeatedly make use of the semi- 
classical limit to explain various signatures of Majorana 
fermions and the topological phase. To demonstrate the 
method, we begin by determining the momentum space distri¬ 
bution of the Majorana fermions in the vortex core and at the 
edge of the system. First, we note that the edge of a topolog¬ 
ically non-trivial superconductor can be seen as an interface 
between the topologically trivial vacuum, and the non-trivial 
bulk. It is therefore a transition from region I to II, similar to 
that along the dashed line in Fig.|^ In the semi-classical limit 
we thus expect the zero-energy spectrum on the edge to be 
related to a gap closing of type D, that is, the Majorana edge 
state will be located at the high symmetry point, k = (tt, tt). 
On the other hand, the vortex core is a region where the most 
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FIG. 2: (Color online). [Top left] Schematic figure of the topologically trivial (I) and non-trivial (II) phases as functions of 1141 and | A |. At the 
topological phase transition following the dashed grey line, the band gap is closed by the formation of a Dirac cone (D). The band gap can also 
be closed by following the full line down to A = 0, at which point the system becomes metallic (M). [I, II, D, M] The two bands closest to the 
Fermi level plotted in the topologically trivial (I) and non-trivial phase (II), at the topological phase transition where a (D)irac cone appears, 
and for A = 0 where the system becomes (M)etallic. 


essential feature is that A —?> 0. The Majorana state in the 
vortex core should therefore be expected to instead be related 
to a gap closing of type M, which produces zero-energy states 
in a circle at a certain radius from k = (tt, tt). 


2. Tight-binding confirmation 

Having deduced the expected behavior from the semi- 
classical limit, we now look at the actual Majorana fermions 
that appear in a fully self-consistent solution. To do so we 
pick the eigenstates corresponding to the two zero-energy lev¬ 
els at point 4 in Fig. [T] and rotate these into the basis where 
the two resulting Majorana fermions and 7 *^^ becomes 
clearly localized in the core and on the edge, respectively. In 
Fig .|^we plot the square of the real space wave function given 
by 


(T 


We also we calculate the momentum space distribution of the 
states 




+ |u 


Mi^2 

k(T 


a 


( 5 ) 


where and are the Fourier transforms of 

and respectively. The result is displayed in the recipro¬ 

cal space plots in Fig. It is clear that the Majorana fermion 
at the edge is mainly built up from k-components at k = 
(tTjTt), while the vortex core Majorana fermion mainly con¬ 
sists of k-components a finite radius away from k = (tt, tt). 


This is in complete agreement with our expectations from the 
semi-classical limit. 

Apart from showing that important features of the Majorana 
fermions can be derived directly from the bulk band struc¬ 
ture by considering the semi-classical limit, this result also 
has important experimental implications; a necessary indica¬ 
tor of the vortex core state being a Majorana fermion is that the 
Fourier transform of the state is distributed along a circle of 
finite radius, and centered at the high symmetry point where 
the bulk band gap closes at the topological phase transition. 
However, we note that this is not an exclusive signature of a 
Majorana fermion, e.g. ordinary Caroli-Matricon-de Gennes 
vortex stateJ^ can be expected to have a similar momentum 
space distribution. The reason being that also these states are 
a consequence of the same local collapse of the superconduct¬ 
ing gap. 


IV. PARTICLE-HOLE SYMMETRY 


Before moving on to further results, we here make a few 
general remarks related to particle-hole symmetry. The BdG 
Hamiltonian for a general and homogeneous bulk supercon¬ 
ductor can be written as a 4 x 4 matrix equation: 


n{k) 


Ho(k) A(k) 
-A*(-k) -H^(-k) 


(6) 


in the basis (ck^, Ck^, c^kf write the eigen¬ 
states as | 7 i,k) = [ Mi/kt From the 

equation T((k) I 7 ) = E\'y) and the complex conjugate of 
I7) = —E I7), it can be seen that these two equa¬ 
tions are in fact the same set of equations. This leads to a 
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Vortex core Edge 
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0 kx 2n 0 kx 2n 


In particular, this implies 

M M M* M* /iA\ 

U-kcr — / ^ ^ k^i^kfT — ko" (I^) 

ly V 

It is therefore clear that the apparent symmetry of the momen¬ 
tum distribution in Fig. is not an accident. Rather, Eq. ( [T0| ) 
guarantees that the momentum distribution of the Majorana 
fermion is inversion symmetric. However, we note that both 
the real space and momentum distribution of the Majorana 
fermions expressed here takes into account both the electron 
and hole components of the eigenstates. This strict symmetry 
is therefore not necessarily present for a physical probe only 
measuring the electron component. 


B. Electronic and BdG expressions 


FIG. 3: (Color online). Majorana fermions in the vortex core and on 
the edge. Real space distribution (top) and Fourier transform of the 
same states (bottom). The vortex core Majorana fermion is mainly 
built up of momentum components at a certain radius away from k = 
(tt, tt), while the Majorana state at the edge is consists predominantly 
of states at k = (tt, tt). 


particle-hole symmetry for eigenstates v and v of opposite 
energy = — i?p. Identification of coefficients in the two 
equations also leads to 

kcr- 

Equation <0 is of particular interest to us because it tells us 
that the states at {E^, k) and {—E^, —k) are related through 

Tpk — ^ ^ ^^pka'Cka' “F — Tp—k' 

a 

This pairwise correspondence is a consequence of the artificial 
doubling of degrees of freedom which occurs when the prob¬ 
lem is treated using a 4 x 4 BdG formulation. The two states 
are therefore not distinct, but rather the occupation of one is 
necessarily accompanied by the deoccupation of the other. 


A. Consequences for Majorana fermions 

The set of provides a complete basis for the solution 
of any not homogeneous problem, with a general eigenstate 
written as 

7^ = XI “Pk7pk = X (^kaCka + , (9) 

i/'k k(T 

where a^k«pk<T and Opk^pkcr- Note 

that the eigenstates are now labeled by their superscript 
A, while the subscript v becomes a summation index. If 
now 7 ^ is a Majorana fermion, it also satisfies the relation 
which using Eqs. (|8]l-(0 leads to ®^7pk = 


The discrepancy just mentioned between what can be stated 
about the Majorana fermion in the BdG formulation and what 
can be seen experimentally reflect a general conceptual diffi¬ 
culty with regard to topological superconductors. The topo¬ 
logical properties are derived from the BdG band structure, 
so the consequences for these properties are most easily un¬ 
derstood in relation to it. However, physical quantities are 
calculated in a way that can obscure their relation to the BdG 
band structure. A common prescription for calculating physi¬ 
cal quantities such as the electronic LDOS, DOS, and spectral 
function involves summations of the form (see e.g. Ref. Il28l ): 

X {\u’;\H{e-e,) + \v’;\H{e + e,)), (ii) 

E^>0 

where p is some index and the summation runs over positive 
energies to avoid over-counting due to the artificial doubling 
of degrees of freedom in the BdG formulation. We note that 
if the 5-function had entered as S{E — E^) in front of the w’s, 
then the u’s would appear as just another orbital, and such 
expressions can be compared straightforwardly with the BdG 
band structure. For this reason it is beneficial to follow an¬ 
other equivalent prescription for calculating physical quanti¬ 
ties. Using Eq. (|7]|, or the same expression transformed to real 
space (depending on the type of index p): it i^ 

possible to replace terms involving u’s in expressions of the 
form in Eq. ( [TT| l, by u terms of the opposite energies. We then 
write 

J 2 {\u';\^ + {\v;\Y)s{e-e,), ( 12 ) 

where the second term in brackets is to be ignored in order for 
the expression to correspond to physical, or electronic, quan¬ 
tities. On the other hand, if the second term is kept, the re¬ 
sulting properties can be straightforwardly related to the BdG 
band structure. The latter expression, which we call BdG- 
type, therefore provides an important conceptual bridge be¬ 
tween the theoretical formulation of (topological) supercon¬ 
ductivity and experimentally measurable quantities. In fact, 
we already used the BdG-type expression in Eq. 0-0 and 
Fig. [3] 
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FIG. 4: (Color online). Eigenstates are classified as edge, vortex 
core, or bulk states according to in which region they are mainly 
located. Eigenfunctions are classified as edge or vortex core states if 
more than 50% of their density are located in the edge or vortex core 
region, respectively. Other states are classified as bulk states. 



V. SPECTRAL FUNCTION 


In Section IIII Al we related the bulk band structure to the 
momentum space structure of the Majorana fermions. Here 
we continue our study of the band structure by solving the real 
space Hamiltonian in Eq. ([^ with a vortex, and then Fourier 
transform the results to arrive at the spectral function. In the 
limit of an infinite homogeneous sample, this is equivalent to 
calculating the analytical bulk band structure. However, when 
the sample is finite and includes defects such as edges and 
vortices, the result will be distorted. The spectral function can 
be calculated as 

ACk,E) = Y, + {\v,ka\'^}) SiE - E,). (13) 


It is useful to divide A(k, E) into parts consisting of edge, 
vortex core, and bulk contributions. For this reason we define 
the state classification functions 


ri 1 ExGX^,a (l^i'xo-P + It'jyxfTp) 

otherwise, 

Cc(^) = 


> 


2 ’ 


Cbii^) = 


1 if J2xeXa,<y (I 

0 otherwise, 

1 ifCe = Ca = 0, 

0 otherwise. 


ttiyxfT I “F It^i/Xfrl ) ^ 2’ 


(14) 


where Xg and Xc are the set of points classified as edge and 
core sites according to Fig. Letting L = e,c,b denote edge, 
core, and bulk, respectively, we define 


ALiK E) = Y, CL{y) (lu.kap + 5{E - E,), 

El,,a 

(15) 


where A(k, E) — 7li(k, E). A peculiarity of this is that 
the numerical zero-energy states, representing the Majorana 
fermions before they are rotated into the Majorana basis, are 
classified as bulk states. This follows since they are located 
both at the edge and in the vortex core, with additional tails 
stretching into the bulk. This conveniently allows us to extract 
from A(k, E) the contribution from the bulk band and Ma¬ 
jorana fermions separately from the contributions from edge 
states and other vortex core states. 

In Fig. 1^ we plot the BdG-type A(k, i?) and Ai(k, i?) 
along ky = TT, at the topologically non-trivial point 4 in Fig. 



Edge Vortex Bulk + MF 
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FIG. 5: (Color online). The BdG-type spectral function A(k, E) 
and its edge, core, and bulk components plotted along ky — tv for 
parameters corresponding to the topologically non-trivial phase at 
point 4 in Fig.[^ A(k, E) is seen to be dominated by contributions 
from the bulk band structure. However, topological edge states, as 
well as vortex core states also contribute. When the edge and or¬ 
dinary vortex core states are subtracted from A(k, E), the resulting 
Ai,(k,E) consists of only bulk bands and the two zero-energy Majo¬ 
rana fermion states. Because the scaling behavior is different for the 
different types of states, the color scale is arbitrary but set to enhance 
the relevant features. 


[2 First of all, we see the two dispersive (at non-zero energy) 
branches of topological edge states in Ae(k, E). Similarly, a 
wealth of vortex core states are visible in Ac(k, E). We note 
that vortex core states predominantly originate from states 
around the band edges. This can be understood by considering 
that at the band edges the DOS is high, which makes it easy 
to hybridize these states into states localized in the vicinity of 
the vortex. We also see that the ordinary vortex core states 
closest to i? = 0 are formed mainly from k-components at a 
distance away from k = (7r,7r), while higher energy vortex 
states are located notably closer to k = (tt, tt). The former, 
which are subgap states, are the ordinary Caroli-Matricon-de 
Gennes states, and their momentum distribution is in agree¬ 


ment with that expected from our discussion in Section III A 


where we mentioned that they, just like the vortex core Ma¬ 
jorana fermions, result from a local (M)-type collapse of the 
bulk gap. In addition, it is clear that Ai,(k,E) contains the 
bulk band structure, as well as faint signals of the Majorana 
fermions at E = 0. 

Although it is possible to image the electronic A(k, E) di¬ 
rectly using e.g. ARPES, it may be difficult to extract the edge 
and vortex core features from such data, due to limited in¬ 
tensity. However, the bulk features such as the Mexican hat 
shaped band structure in Fig.j^scale with the area of the bulk, 
and can therefore be imaged directly with ARPES. On the 
other hand, it should be feasible to image the vortex core con¬ 
tribution by Fourier transforming data from local probes such 
as STS. Note, however, that it is not the BdG-type spectral 
function in Fig.|^that is physically measured. Rather it is the 
corresponding electronic spectrum which results from drop¬ 
ping the hole-part in Eq. ([T5]|. Having seen how various fea¬ 
tures of the spectral function can be understood by splitting it 
decomposing into edge, vortex core, and bulk states contribu- 
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tions, we define Ah_|_c(k, E) = Ai,{k, E) + Ac(k, E). This is 
the spectral function with contributions from edge states ex¬ 
cluded, which are artificial effects introduced by the finites 
size of our sample. We will use this spectral function in the 
rest of this article. 


VI. LOCAL DENSITY OF STATES 

Having investigated signatures in momentum space, we 
now turn to a discussion of the LDOS, given by 

= + (16) 
1^,(7 

We will see that a comparison between the LDOS and spec¬ 
tral function is helpful for revealing important signatures of 
the topological phase. In addition we also compare with the 
bulk DOS calculated from a pure bulk solution with an equiv¬ 
alent I A|. In Fig. I^the LDOS, bulk DOS, and Ab+c(k, E) are 
plotted side by side at the points labelled 1-4 in Fig. [T] Both 
the electronic and BdG-type results are shown here. In what 
follows we use the BdG-type results for the interpretation, but 
the electronic results are to be considered when comparing 
to experiments. When the Zeeman field is zero (1), the band 
structure can be seen to consist of two Rashba spin-orbit split 
parabolas crossing at a high symmetry point. The band edges 
shows up as sharp features in both the LDOS and DOS. Vor¬ 
tex core states below the band edges can also be observed. As 
the Zeeman field is turned on (2), the two bands are split off 
from each other, and the upper band becomes irrelevant for 
the low-energy spectrum making the system effectively ’’spin¬ 
less”. Still vortex core states form below both band edges. 
Inside the F region (3), a vortex core state appears at or close 
to E = 0, and there gives rise to a non-Majorana fermion 
zero bias peak. It is clear from the spectral function that the 
topologically non-trivial phase has not yet been entered, as 
the two parabolic bands closest to the Fermi level do not give 
a Mexican hat shaped spectrum. However, this can not be 
seen from the LDOS at = 0 alone, and care thus has to 
be taken to not casually interpret every zero-energy peak in 
the LDOS spectrum as a Majorana fermion. We further note 
that all vortex core states in these first three cases are asso¬ 
ciated with faint signals in the spectral function around the 
high-symmetry point and below the band edges. This is par¬ 
ticularly important for the zero-energy states in region F, as 
this make their momentum distribution clearly distinct from 
those of the Majorana vortex core states. 

Once the two low-lying bands becomes inverted and ac¬ 
quires Mexican hat like forms, the topologically non-trivial 
phase is entered, as exemplified by point (4). Because of the 
comparatively rather flat nature of the Mexican hat like struc¬ 
ture, the band edges acquires a large density of states, which 
generates a strong intensity in the (L)DOS. The observation of 
strong band edges in the LDOS spectrum therefore provides 
evidence for the system being in the topologically non-trivial 
phase. However, this only distinguishes the topologically non¬ 
trivial phase from the trivial phase, and not from a conven¬ 
tional superconductor where strong band edges also appears. 


We also note that the bulk DOS clearly shows that the band 
edge in fact consists of two edges, one from the Mexican hats 
lower edge away from k = (tt, tt), and one from its edge at 
higher energy at k = (tt, tt). This provide a clear distinction 
also from the band edge behavior of a conventional supercon¬ 
ductor. Additionally, we see that the zero-energy Majorana 
fermion is accompanied by a wealth of other vortex states, 
showing up as a distinct x-shape structure in the low-energy 
LDOS. These subgap states are the ordinary Caroli-Matricon- 
de Gennes states already discussed. 

Another important piece of evidence for the topologically 
non-trivial phase is provided by an apparently larger band gap 
inside the vortex core. From the LDOS it appears as if the 
band gap is increased around the vortex, leaving a strong vor¬ 
tex core signal above the bulk band edge. This is due to the 
fact that the vortex core states now behave differently than the 
vortex core states observed so far, which always appear below 
the bulk band edge. We can understand this as a consequence 
of the two different band edges of the Mexican hat shaped 
bulk band. The lower band edge away from k = (tt, tt) col¬ 
lapses in the core because it is due to superconductivity. How¬ 
ever, the second band edge does not collapse. Rather it can be 
seen from Eq. © that at k = (tt, tt) the energy is given by 

En = ± ^|I41 — ■ It further clear that 

in the topologically non-trivial phase the square root term is 
necessarily smaller than |I4|. Thus, as |A| —)• 0 inside the 
vortex core, the second band edge is pushed away from the 
Fermi level. This result in a set of vortex core states above 
the band edge. For more details of this process we provide a 
zoomed-in plot of the LDOS at point (4) in Fig. 13 First of all, 
the Majorana fermion is clearly visible as a state atE = 0, and 
has a clear energy-separating from the other Caroli-Matricon- 
de Gennes states, making it possible to resolve. Shown is also 
the BdG-type DOS, obtained by summing the bulk and vor¬ 
tex core DOS of the self-consistent solution, the vortex core 
DOS for the self-consistent solution, as well as the DOS of the 
pure bulk solution. A comparison between the self-consistent 
bulk-Hcore DOS and the DOS of the pure bulk solution clearly 
show how the lower band edge essentially collapses around 
the vortex and gives rise to sub gap states. On the other hand, 
the upper band edge is pushed up in energy, giving rise to the 
vortex core states above the band edge, most clearly visible as 
a half x-shape structure above the band edge in the electronic 
spectrum. 

We finally also point out that the Majorana fermion is well 
localized in the center of the vortex core, and that the arms of 
the x-shape structure formed by the subgap vortex core states 
meet at the center of the vortex. We put this in contrast to re¬ 
cent STS experiments, where a spatially extended zero-energy 
peak, together with the arms of the x-shape structure pointing 
towards zero a finite distance away from the center of the vor¬ 
tex, were assumed to provide evidence of Majorana fermions 
in a superconductor-topological insulator heterostructure.l^ 
Our results clearly show that Majorana fermions do not pro¬ 
duce such signatures in a spin-orbit coupled semiconductor. 
However, it should be noted that these experiments were done 
on an s-wave superconductor coated by a topological insula¬ 
tor, and may therefore be sufficiently different that a direct 
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FIG. 6: (Color online). LDOS for a cut through the vortex core (left), bulk DOS (center), and spectral function at ky = 7r(right), for the four 
sample points shown in Fig. [T] using both the electronic and BdG-type expressions. The position of the vortex core states in the LDOS and 
spectral function are indicated in the DOS plots by dotted lines, and appears in grey if they are associated with the two higher-energy bands. 


Electronic BdG 



FIG. 7: (Color online). Same LDOS as in (4) in Fig. but with 
higher resolution around E = 0. The vortex core states are now 
seen to be clearly separated by finite energies. The DOS is plotted 
to the right. Leftmost DOS consists of bulk and core contributions 
from vortex calculation. The middle DOS contains only vortex core 
states. Rightmost DOS is from a pure bulk calculation. 


comparison is not possible. 


VII. PAIR AMPLITUDES 

It has recently been suggested that the appearance of Ma- 
jorana fermions i s clo sely related to the presence of odd- 
frequency pairing.1^21111 in addition, the appearance of uncon¬ 
ventional pair a^litudes in vortices has attracted great in¬ 
terest in general.®^ For this reason we here also provide a 
detailed investigation of the superconducting pair amplitude 

i^(R,r, (T, cr') =(cR4.r,crCR_r,a') 

— 'y ' '*^i/,R+r,tT^i',R.—) ( 1 '^) 

which we decompose into singlet and triplet components 
F*(R, r) (F(R, r, f, i) - F{R, r, i, f)) , 

F™=i(R,r) =F(R,r,t,t), 

(R, r) = i (F(R, r, t4) + r, t)), 

F™=-i(R,r)=F(R,r,i,i). (18) 

Here R is the center-of-mass coordinate for a pair of electrons, 
while r and — r points to the two electrons. The pair ampli¬ 
tude can further be classified as s- p- d-wave, and so forth, ac¬ 
cording to its angular dependence around the center-of-mass 
coordinate. For this purpose the pair amplitude is further pro¬ 
jected onto e**® to obtain the s-, p+-, p_-, d+-, and d_-wave 
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FIG. 8: Orbital wave functions defined by the projection operators 
in Eq. CD- The values on the sites are the complex conjugate of 
the value of the projector when r points from the central site to the 
corresponding site. 

pair amplitudes for I = 0,1, —1, 2, —2, respectively. In gen¬ 
eral, there is a certain degree of freedom in how to perform 
this projection, because the pair amplitude also has a radial 
dependence. In practice, however, the pair amplitude is only 
expected to be sizable within a few coherence lengths, and it 
is therefore only the components of the pair amplitude which 
corresponds to small r that is of interest. Focusing on the 
pair amplitude components for minimal r (without becoming 
trivially zero) for the various orbital moments, the appropriate 
projectors are 

^^a., 

^P± ~ ^ , 

a^A 

=5 (19) 

° bGB 

with A = {(1,0), (0,1), (-1,0), (0,-1)}, i? = A u 
1(1,1), (—1,1), (—1, —1), (1, —1)}, and 0(r) the angular co¬ 
ordinate of r. The orbital wave functions onto which the pair 
amplitude is projected are displayed in Fig. The projected 
pair amplitudes can now be written as 

FS(R)^Y.PoF^{R,r), (20) 

r 

where we study O G {s, Sext,P±, d±} and S G (s, 1,0,-1}, 
while r runs over all possible lattice vectors. 

A. Even-frequency pairing 

In Fig. we plot both the maximum value as a function of 
the Zeeman field, as well as representative real space profiles 
of the pair amplitudes. All pair amplitudes are normalized by 
the bulk value of the conventional s-wave (singlet) pair am¬ 
plitude at zero Zeeman field. The conventional s-wave pair 
amplitude has been excluded from the plots as it is directly 
given by the relation |A(R)| = | — 14c.P(R-, 0)t)l- From 
the real space profiles it is clear that three types of behavior 
can be identified. Namely, the pair amplitudes can be seen as 
originating in the bulk, vortex core, or region surrounding the 
vortex core, which here will be referred to as the pre-core re¬ 
gion. All pair amplitudes have their maximum value in one of 


these three regions, and decays to zero in the other two. It is 
also clear that the three most important pair amplitudes in the 
bulk, apart from the s-wave (singlet) pair amplitude, are the 
extended s-wave (singlet), (m = — 1), and p_ (m = 1) 
components. These all have a total angular momentum z-axis 
project Jz = 0, as expected if angular momentum is to be 
conserved. The appearance of extended s-wave (singlet) is 
a direct consequence of conventional s-wave superconductiv¬ 
ity; it only gives further details about the size of the pair am¬ 
plitude at a finite radius r. The p-wave pair amplitudes in the 
bulk appear because of the Rashba spin-orbit interaction and 
is in agreement with previous results.l^ Our results do how¬ 
ever provide further details showing that the p+ (m = — 1) 
component is largest when the Zeeman field is positive, while 
p_ (m = 1) dominates for a negative Zeeman field. This can 
be understood since essentially all electron levels are occupied 
in a lightly hole-doped semiconductor. Thus the band cross¬ 
ing the Fermi level is the band which is pushed up in energy 
by the Zeeman field and contains predominantly spins anti¬ 
aligned with the Zeeman field. This is also the band that is 
gapped by superconductivity, explaining the overweight of su¬ 
perconducting pairing among spins anti-aligned with the Zee- 
man field. We note that the argument would be reversed for a 
lightly electron-doped semiconductor. 

Next, we note that the existence of two p_ states, in the 
vortex core (m = 0) and in the pre-core region (m = — 1), 
are strongly correlated with the onset of the wide vortex core 
region F. For negative 14, for which no vortex core widening 
occurs, we already have p_ (m — 1) preferred in the bulk. 
The rotation direction of the p_ orbital part agrees with the 
rotation direction of the vortex. For positive Vz, on the other 
hand, the (m = —1) state is preferred in the bulk but it 
has an orbital motion directed opposite to that of the vortex. 
This appears to lead to a widening of the vortex core. Finally, 
we note that in the bulk and vortex core the total angular mo¬ 
mentum is Jz = 0 and Jz = —1, respectively. The angular 
momentum Jz = —1 of the vortex core pairing can be ex¬ 
plained by the vortex winding n = — 1 being absorbed into 
the orbital part of the pair amplitude.HSl On the other hand, 
the pre-core pair amplitudes respect neither the bulk nor the 
vortex core angular momentum conservation rule. 


B. Odd-frequency pairing 

Having described the behavior of the regular pair ampli¬ 
tudes we now turn to an investigation of odd-frequency pair 
amplitudes. The possibility of odd-frequency pairing arise 
when r is considered to not only be a coordinate in space, but 
also in time, which we denote by r = (r, f/2). The ordinary 
requirement on the pair amplitude to be odd under the simul¬ 
taneous interchange of position and spin (cR+f_crCR-f,cr') = 
~ (cR-f ,(t' CR+f ,dr) is then transformed to the requirement that 
the pair amplitude is odd under the simultaneous interchange 
of position, spin, and time (frequency). The pair amplitudes 
investigated so far are all even in frequency, as only the even 
frequency components can be non-zero for f = 0. To identify 
pair amplitudes which are odd in time, we have to consider 
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FIG. 9: (Color online). Even frequency pair amplitudes. [Left] maximum absolute value of the pair amplitudes as a function of 14. [Right] 
representative profiles for of the pair amplitudes along the a;-axis, through the vortex core. All profiles are plotted for 14 = 0.41, just to the 
right of the transition from I to F, except the m = 1 profiles where 14 = —0.5, but all pair amplitudes retain their profile characteristics as 14 
is varied. All values are normalized by the bulk value of the spin-singlet s-wave pair amplitude at zero Zeeman field, Fg{bulk, 0). 


the derivative of the pair amplitude with respect to timd^Ull 


dF(R, r, CT, cr') d 
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',R+r,fT^R-— 
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Preforming the decomposition of the total pair amplitude into 
its spin and orbital components as outlined in Eq. ( fTTpOl l al¬ 
lows us to also consider the amplitudes dFQ (R, t) / dt. 

In Fig. 10 the odd frequency equivalent of Fig. is dis¬ 
played. Most notable is the appearance of an s-wave (m = 0) 
pair amplitude in the bulk when the Zeeman field is non-zero. 
Likewise, the extended s-wave (m = 0) component is found 
to have a similar behavior. All the pair amplitudes associated 
with the bulk have = 0. 

Next, we note four non-zero pair amplitudes in the core. 
These are the s (m = —1), extended s {m = —1), p- (sin¬ 
glet), and d- (m = 1), which all have = —1. The exis¬ 
tence of these can once again be understood as a consequence 
of the vortex winding being rotated into the pair orbital part. 
However, although small, an anomalous d+ {m = 1) compo¬ 
nent also appear in the core, violating the otherwise seemingly 
perfect agreement with vortex core pair amplitudes having 
Jj = — 1. This is due to the fact that each site has four nearest 
neighbors, while four of the sites included in the projection 
onto the d-wave pair functions are second-nearest neighbors. 
On the nearest neighbor sites so with re¬ 

spect to these site alone the angular momentum can only be 
defined modulo 4, and consequently — 3 mod 4 = — 1. 
This is an example of how a radial variation in the pair am¬ 
plitude, in this case the d- (m = 1) component, can give rise 
to seemingly angular momentum violating components when 
calculated on a discrete lattice. Unfortunately it is impossible 


to define a completely satisfactory projection procedure on a 
discrete lattice. Finally, we note the existence of the pre-core 
type pair amplitudes s (m = 1), extended s (m = 1), p+ 
(singlet), and d+ (m = —1). These, just as the pre-core type 
even-frequency p±, all have satisfying neither = 0 nor 
Jz = —1. They do, however, all have in common that they 
either have a bulk or core counter-part differing only in spin 
or orbital rotation direction. We therefore interpret these as 
being secondary in nature, being induced from their bulk and 
core counter-parts. 

Before ending the discussion of pair amplitudes we also 
note that although the odd frequency s-wave (m = —1), as 
well as even frequency (m = —1) and p_ (m = 1) pair 
amplitudes all increase as function of the Zeeman field to¬ 
wards the non-trivial phase, they contain no specific signature 
of the non-trivial phases themselves. We put this in contrast to 
recent results relating the existence of odd-frequency s-wave 
pair amplitudes with the appearance of Majorana fermions in 
ID.I^ Likewise, the various pair amplitudes which appear in 
the vortex core provides little clue to the topological phase or 
existence of Majorana fermions, but appear in the whole phase 
diagram. On the other hand, several of the pair amplitudes 
changes abruptly at the Jo/' transition, thereby reflecting 
the sudden widening of the vortex core. A summary of the 
described pair amplitudes is provided in Table [I] 


VIII. SUMMARY 


We have in this work investigated the local density of 
states (LDOS), band structure, and superconducting pair am¬ 
plitude for signatures of the non-trivial topological phase 
and Majorana fermions in vortex cores in spin-orbit coupled 
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FIG. 10: (Color online). Same as Fig.|^ but for odd frequency. All real-space profiles are for 14 = 0.38. 


Orbital 

Frequency 

Spin 

J. 

Origin 

Max (%) 

P+ 

even 

m = 1 

2 

pre-core 

2 

P+ 

even 

m = 0 

1 

pre-core 

1 

^ext 

even 

singlet 

0 

bulk 

10 

p+ 

even 

m = —1 

0 

bulk 

14 

p- 

even 

m = 1 

0 

bulk 

14 

p- 

even 

m = 0 

-1 

core 

7 

p- 

even 

m = —1 

-2 

pre-core 

5 

d- 

even 

singlet 

-2 

pre-core 

1 

d+ 

odd 

m = 1 

3 

core 

2 

s 

odd 

m = 1 

1 

pre-core 

1 

^ext 

odd 

m = 1 

1 

pre-core 

1 

p+ 

odd 

singlet 

1 

pre-core 

2 

d+ 

odd 

m = —1 

1 

pre-core 

1 

s 

odd 

m = 0 

0 

bulk 

41 

^ext 

odd 

m = 0 

0 

bulk 

5 

S 

odd 

m = —1 

-1 

core 

6 

^ext 

odd 

m = —1 

-1 

(pre-)core 

6 

p- 

odd 

singlet 

-1 

core 

10 

d- 

odd 

m = 1 

-1 

core 

5 


TABLE I: Even and odd pair amplitudes, ordered according to their 
Jz quantum number. Jz = 0 and Jz = — 1 pairing is seen to orig¬ 
inate in the bulk and core, respectively, while other pair amplitudes 
only appears in the pre-core region. The d+, odd frequency, m = 1 
pair amplitude is anomalous in that it originates in the core region 
although it has Jz = 3. As explained in the text, it is an artifact of 
the pair amplitudes being calculated on a lattice. Likewise, the Ssxt, 
odd frequency, m = — 1 pair amplitude is anomalous in the wide 
vortex core region, as it there originates in the pre-core region rather 
than in the core, in spite of having Jz = —1. The maximum relative 
strength as compared to the s-wave, even frequency, spin-singlet pair 
amplitude at zero Zeeman field is tabulated at the right. 


semiconductor-superconductor heterojunctions. A necessary 
indicator of a zero-energy vortex core state being a Majo- 
rana fermion has been identified to be a momentum distribu¬ 
tion centered at a finite radius away from the high symme¬ 
try point k = (tt, tt), [k = (0, 0) in the case of a lightly 
electron-doped semiconductor]. Moreover, the vortex Ma- 
jorana fermion and finite-energy Caroli-Matricon-de Gennes 
vortex states are found to be well separated in energy and in 
the topological phase they together form a characteristic x- 
shape structure in the subgap LDOS when scanning through 
the vortex core. The Majorana mode is very well-localized in 
the center of the core, while the finite-energy states disperse 
further out from the center, although the x-shape structure is 
still centered at the core center. 

Furthermore, we show that a clear signature in the spec¬ 
tral function of the topological phase itself is the Mexican 
hat shaped band structure, which also gives rise to double 
band edges, very clearly visible in the DOS due to their high 
concentration of states. These double band edges also give 
rise to the existence of a second class of vortex core states, 
distinct from the ordinary Caroli-Matricon-de Gennes vortex 
core states and the Majorana fermion. These vortex states ap¬ 
pear beyond the superconducting gap and forces locally larger 
band gap in the vortex core region, both features that are ex¬ 
perimentally measurable. 

Finally, we have also investigated the superconducting pair 
amplitude, showing that multiple pair amplitudes with total 
Jz = 0, both even and odd in frequency, develops in the bulk 
because of the finite spin-orbit interaction and magnetic field. 
In the vortex core we instead find pair amplitudes which have 
a total Jz = — 1, where the vortex core momentum (— 1) has 
been rotated into the orbital part of the pair amplitudes. De¬ 
spite multiple unconventional pairing amplitudes developing 
in the core, we find no amplitude that signals the onset of non¬ 
trivial topological order. Specifically, the appearance of a Ma¬ 
jorana fermion does not imply any noticeable increase in the 
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odd-frequency pairing. However, we find a strong correlation 
between the transition from narrow to a wide vortex core and 
the development of even-frequency p_ components. In sum¬ 
mary, these results provide multiple specific characteristics 
for the non-trivial topological phase and its vortex Majorana 
fermion in spin-orbit coupled semiconductor-superconductor 
heterostructures. These distinct indicators provide both added 
physical understanding of the topological phenomena in these 
heterostructures and are in many cases directly experimentally 
measurable. 
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